Many-body Landau-Zener dynamics in coupled ID Bose liquids 
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The Landau-Zener model of a quantum mechanical two-level system driven with a hnearly time 
dependent detuning has served over decades as a textbook paradigm of quantum dynamics. In 
their seminal work [L. D. Landau, Physik. Z. Sowjet. 2, 46 (1932); C. Zener, Proc. Royal Soc. 
London 137, 696 (1932)], Landau and Zener derived a non-perturbative prediction for the transition 
probability between two states, which often serves as a reference point for the analysis of more 
complex systems. A particularly intriguing question is whether that framework can be extended to 
describe many-body quantum dynamics. Here we report an experimental and theoretical study of a 
system of ultracold atoms, offering a direct many-body generalization of the Landau-Zener problem. 
In a system of pairwise tunnel-coupled ID Bose liquids we show how tuning the correlations of the 
ID gases, the tunnel coupling between the tubes and the inter-tube interactions strongly modify 
the original Landau-Zener picture. The results are explained using a mean-field description of the 
inter-tube condensate wave-function, coupled to the low-energy phonons of the ID Bose liquid. 



Systems of ultracold bosonic atoms in optical lattices 
provide clean and highly tunable laboratories for inves- 
tigating strongly correlated quantum states in low di- 
mensions. Prominent examples in one dimension are the 
conversion of weakly interacting superfluids into strongly 
correlated Tonks-Girardeau gases [1-4] or Mott insulat- 
ing states [5] including the measurement of hallmark cor- 
relations [6, 7] and transport properties [8]. 

Next to the investigation of equilibrium phases, current 
ultracold atomic systems offer an entirely new perspec- 
tive on low-dimensional quantum liquids, allowing one 
to focus on their non-equilibrium dynamical behaviour. 
Compared to the equilibrium properties, the theoretical 
understanding of dynamics in such systems is far more 
rudimentary and is currently a topic of intense study [9- 
15]. A few experiments utilized sudden changes in system 
parameters in order to induce a quantum quench. This 
allowed to study fundamental questions including ther- 
malization in low-dimensional systems [16] and interfer- 
ometry using many-body dynamics [17-19]. 

The limit of adiabatic dynamics, however, is equally in- 
teresting and offers a more likely opportunity to observe 
universal physics [20-22]. The basic prototype of such 
phenomena is the Landau-Zener (LZ) problem, which 
describes a sweep through an anti-crossing of two lev- 
els [23, 24]. The two states can be thought of as those 
of a single particle in two potential wells, which differ 
in energy by a time-dependent detuning A = a • t and 
are tunnel-coupled with strength J. The celebrated LZ 
formula 

F,z(J,a) = l-e-2-^ (1) 
pertains to the probability for a particle initiated in the 
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bottom well at i = — oo to end up in the opposite well at 
t — +00, following a linear sweep of the detuning at the 
rate a. The solution of the LZ problem often serves as 
a starting point for the analysis of more complex prob- 
lems ranging from sweeps across a Feshbach resonance 
along the BEC-BCS crossover [25-27] to transport and 
dissipation in mesoscopic systems [28, 29]. 

Here we report on an experimental and theoretical in- 
vestigation of the controlled LZ dynamics in a strongly 
interacting many-body setting of two tunnel-coupled ID 
quantum liquids. Next to a strong dependence of the 
LZ transition probability on the inter-particle interac- 
tions and the condensate fraction, we find a complete 
breakdown of adiabaticity in the system for the case of 
a LZ sweep in the excited state, where the transfer fi- 
delity is decreased rather than increased at low sweep 
rates a. We attribute this breakdown to self-trapping of 
the meta-stable condensate near zero detuning and sub- 
sequent decay by phonon emission into the Bose liquid. 
Strong correlations, which develop in the quantum liquid 
upon approaching the Mott-insulating phase arc found 
to suppress the self-trapping and decay mechanism, thus 
restoring the adiabatic transfer. 

Our system realizes the controlled LZ dynamics in an 
array of pairwise coupled potential tubes oriented along 
the z-axis, which are created using a two-dimensional op- 
tical lattice potential with a supcrlattice potential along 
the orthogonal x-direction (see Fig. la). The latter is 
formed by two coUinear retro-reflected lattice standing 
waves with a periodicity ratio of 2 and tunable relative 
phase [30, 31]. All lattice depths are given in units of the 
respective recoil energies E^. = h? / [2mB.b^],) with rriRb 
being the atomic mass (see Methods). Our experiments 
begin by loading a BEG of about 9 x 10^ ®''Rb atoms 
in the \F = l^mp = —1) state from a magnetic trap 
into the 2D optical lattice, translating to about 100 par- 
ticles in the central tubes (see Methods). At the outset 
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FIG. 1; Landau-Zener sweeps in a quantum ladder. a 
Schematic sketch of the experimental setup. Pairs of ID Bose 
liquids are tunnel coupled through a potential barrier along 
the a;-direction. In this system, the potential offset (bias) A 
between the tubes, the inter- and intra-tube kinetic energies 
J and J|| , the chemical potential ^ and the correlations prop- 
erties of the ID gases characterized by the Luttinger param- 
eter K can be controlled, b Measured relative populations 
in right tube n_R after ground state sweeps as a function of 
the inverse sweep rate /a with single atoms (black cir- 
cles), single atoms and pairs (grey circles) and with ID Bose 
liquids (open circles). The insets show typical band mapping 
images from which the populations in the left and right well 
are deduced. 



only one tube of each coupled pah', the left tube, is filled 
with atoms. The remaining empty right tube is far de- 
tuned with initial detuning A^, either below (A^ > 0) or 
above (Aj < 0) the left tube. Then the detuning A is 
changed linearly in time by changing the relative phase 
of the superlattice, causing particle transfer via tunnel- 
ing between the tubes. After the sweep, we measure the 
relative population {nji) in the left (right) tubes re- 
spectively, by applying a site-resolving mapping method 
[31, 32] (see insets Fig. lb). 

The experimental setup allows us to study how the 
transfer rate between the tubes depends on the many- 
body correlations of the ID Bose gases, and on the ratio 
between particle interactions and the inter-tube tunnel 
coupling. The former can be controlled via an additional 
lattice along the tube axis (the axial ^-lattice) , while the 
latter can be adjusted via the potential-barrier height 
set by the depth of short-wavelength lattice along the 
transverse x-direction. 



I. GROUND-STATE SWEEP 

We begin our discussion with the case of ground-state 
sweeps, where A^ < such that the filled left tube is 
the one with lower energy. The comparison of transfer 
efficiencies {nn) for pairwise coupled dot-like lattice sites 
and one-dimensional tubes as a function of sweep rate 
is shown in Fig. lb. The former scenario is realized 
by adding a deep lattice along the axial z-direction. To 
quantify the sweep fidelity with a single number, we de- 
fine the characteristic rate ac as the sweep rate to achieve 
a transfer efficiency of 1 — e^^ « 63% with respect to 
the maximal measured value. Therefore, larger ac cor- 
respond to higher transfer rates. In the case of single 
atoms in the dot-like sites, the result almost matches 
with the LZ formula Eq. (1). Here, ac = 27rJ^ sets the 
natural scale for the sweep rate. If the filling is increased, 
we find an enhanced value of a^ This enhancement is 
even more pronounced in the case of the pairwise coupled 
tubes filled with up to 100 atoms (see Fig. lb). 

The increase of transfer rate with atom number can al- 
ready be understood within the simple two-mode model 
of a "double dot" occupied by strongly interacting par- 
ticles [33, 34]. For the case of strong interactions within 
a mode compared to the tunnel coupling, the LZ tran- 
sition is split into N avoided level crossings, each corre- 
sponding to the transfer of one particle with a coupling 
enhanced by Bose statistics (see Supplementary Informa- 
tion). Note that interactions are crucial to observe this 
enhancement: a non-interacting BEC would result in the 
same transfer rate as is the case for a single particle. 

From now on we focus on the the case of large parti- 
cle numbers. Then an effective ID Bose liquid forms in 
the filled tubes, and we study how the transfer efficiency 
to the empty tubes depends on the properties of that 
liquid. We control the nature of the many-body system 
by tuning the strength of both the transverse z-lattice 
and axial z-lattice potentials. Both parameters strongly 
influence the characteristic sweep rate ac in the LZ prob- 
lem (see Fig. 2): ac is found to increase with increasing 
strength of the transverse x-lattice potential while it de- 
creases when increasing the axial ^-lattice depth. For 
very large z-lattice depth, the system is eventually driven 
across the Mott transition in the ladder system [35] and 
deep in the Mott-insulating regime for unity filling, the 
single-particle result is restored. 

In order to explain this behaviour, we consider an in- 
finite homogeneous Bose-Hubbard ladder with an aver- 
age of one particle per rung (see Methods). We treat 
the inter-tube coherence by a mean-field approximation. 
That is, we take the quasi-condensate in the ladder to 
be in a fixed linear combination of left and right tube 

= V'L^L + V'fl^lii- This leads to a reduced single-chain 
Hamiltonian If^ff [^l, V'i?,], which depends on the varia- 
tional parameters i/'l and ipR (Eq. (6) in Methods). Min- 
imization of the ground-state energy E(ip[^,ipf{), subject 
to the constraint [V'lI^ + \4'r\'^ = 1 defines a non-linear 
LZ problem, that is analogous to the two-mode non-linear 
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FIG. 2: Influence of interactions and correlations on the char- 
acteristic sweep rate. Experimentally measured scaled char- 
acteristic sweep rate ad as a function of the trans- 
verse (a) and axial lattice depth (b). The insets (c,d) show 
the theoretical predictions using mean-field theory (Eq. (2)) 
as described in the text. 



pling strength -y/gJ. Here ac = 27rqJ'^, with q given in 
ref. [41] and in the Supplementary Information as a func- 
tion of the effective non-linearity rj. 

The characteristic rates ac calculated in this way are 
shown as insets in Fig. 2. Their dependence on the 
transverse x-lattice and the axial z-lattice depth accounts 
well for the trends observed in the experiments. Both of 
these trends stem from the change of the effective non- 
linearity rj. An increase in the potential barrier between 
the tubes (z-lattice depth) increases rj and reduces the 
curvature of the ground-state energy in its dependence 
on detuning, thus allowing for faster sweeps [36] (corre- 
sponding ultimately to the split resonance case discussed 
above for rj :s> 1). Increasing the z-lattice depth along 
the tubes first increases 77 via an increase in the effec- 
tive U. For even larger lattice depths, however, the con- 
densate is significantly depleted when driven towards the 
Mott transition of the quantum ladder. Deep in the Mott 
regime all Bose enhancement in the tunneling process is 
lost and one essentially recovers the single-particle LZ 
results. The fact that the calculated ac goes to 1 imme- 
diately at the transition is an artifact of the mean-field 
theory (Eq. (2)). Quantum fluctuations beyond mean- 
field theory give rise to on-site particle number fluctu- 
ations and coupling between the sites even in the Mott 
phase [? ]. This will smoothen out the change of ac 
across the transition, saturating to the asymptotic value 
ac/27rJ^ — 1 deep in the Mott-insulating regime. We 
attribute the slightly higher saturation value seen in the 
experiment (ac ~ 1.5) to the fraction of doubly occupied 
sites that are present in the system. 



Schodinger equation in ref. [41, 42]. 

To solve for £'(^/;l, i/'ij) we employ a mean-field theory 
of the single-chain Hamiltonian (see Eq. (6) in Methods), 
which leads to decoupled local Hamiltonians on rungs of 
the ladder: 

^^MF = -J||(Ml + $*a,) + t/(|^L|^ + |^/?|')aI4a.a. 
+ [A(|^i|2 - |^^|2) _ J(^2^^ + c.c.)] a|a,. (2) 

Here $ = {cii) is the mean-field condensate. While 
this approximation does not capture the subtle long- 
range correlations of the one-dimensional system, it is 
rather good for the energetics. In the weak coupling 
limit (fli) ~ -y/n, and Eq. (2) reduces to the usual two- 
site Gross-Pitaevskii functional [42]. At stronger inter- 
actions, the condensate $ is depleted until it vanishes at 
the transition to the Mott-insulating state in the ladder. 

Within the mean-field Eq. (2), the time dependence of 
the two modes ipL and ipji is governed by a non-linear 
Schrodinger equation similar to ref. [41, 42], with the 
effective non-linearity rj « U^'^/J (see Supplementary 
Information). We can therefore apply an adiabatic per- 
turbation theory [41], in order to obtain a modified LZ 
formula P/j = 1 — exp{—2TTqJ^ /a) with an effective cou- 



II. INVERSE SWEEP 

We now turn to the investigation of the inverse LZ 
sweep. Namely, we initialize the system with > 
such that the filled left tube is the one with higher energy. 
We follow the population of the tubes as the detuning is 
gradually changed until the relative potential between 
the tubes is reversed. 

In Fig. 3a we plot n/j as a function of the final detuning 
Ay for a relatively slow sweep with 27rJ^/|Q;| = 2.1(1). 
While the ground-state sweep shows a continuous trans- 
fer of population, the inverse sweep exhibits a slightly 
reduced transfer for A/ > and an abrupt breakdown 
for A/ < once the resonance is crossed. From then on, 
particles are transferred back to the left tube which is 
now lower in energy. We will argue below that the dra- 
matic difference between the two cases stems from the 
highly non-equilibrium nature of the inverse sweep com- 
pared to the near-equilibrium conditions of the ground- 
state sweep. In Fig. 3b we illustrate the two mechanisms 
responsible for the breakdown of adiabaticity: 1) axial 
phonons which provide a relaxation mechanism and 2) a 
loop structure of the mean-field condensate energy, which 
leads to self-trapping. 

Fig. 4a shows a more detailed study of the adiabatic- 
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FIG. 3: Adiabaticity breakdown in the inverse LZ sweep, a 
Comparison of the transfer efficiency ur from tire filled to 
the initially empty tube as a function of the final detun- 
ing for a ground state sweep (black dots) and an inverse 
sweep (red dots). The sweeps were carried out at a con- 
stant rate 2i\j'^ /\a\ = 2.1(1) and a;-lattice depth of 15 . 
For the inverse sweep a sudden breakdown in transfer fi- 
delity is observed beyond resonance, b Adiabatic energy lev- 
els for ground state sweep (black curve) and inverse sweep 
(red curve). For strong repulsive interactions a loop develops 
in the condensate energy level structure [42]. In an inverse 
sweep the condensate becomes self-trapped as it follows the 
path along the upper branch until reaching the spinodal point 
in the outer edge of the loop (red dot). Beyond this point it 
becomes unstable and decays to the lower branch, while emit- 
ting phonon excitations along its axial direction. The dashed 
blue curve denotes the energy levels of the non-interacting 
problem. 



ity breakdown for different axial z-lattice depths. We see 
that at fast sweeps (27rJ^/|a| < 1), hr increases with 
decreasing sweep rate as expected from the standard LZ 
problem. In this regime, the curves corresponding to dif- 
ferent z-lattice depths follow the typical universal scal- 
ing with 2'iTj'^/\a\. By contrast, in the slow sweep regime 
27r J^/laj > 1, this scaling is no longer obeyed. Moreover, 
we find the transfer efficiency n/j to actually decrease 
with decreasing sweep rate for sufficiently low z-lattice 
depths, indicating a complete breakdown of adiabaticity. 

Such a decay can be understood as follows. Particles 
occupying the upper state close to A = 0, i.e. the high- 
energy superposition of left and right tubes, form a filled 
ID Bose liquid. The orthogonal low-energy superposition 
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FIG. 4: Influence of intra-tube correlations on the adiabaticity 
breakdown, a Measured transfer efficiency na as a function 
of the scaled inverse sweep rate for different axial z-lattice 
depths. The x-lattice depth was 12 E^, and the z-lattice 
depths were set to S^, 2 i;^, 4 E^, 7 S^, 10 and 20 
Er, respectively. The inset shows the transfer efficiency ex- 
tended to slower sweeps, b Transfer probability as a function 
of scaled inverse sweep rate computed using an effective Lut- 
tinger liquid theory (see Eqs. (3,4)) valid at weak non-linearity 
(see text). 



state forms an effective empty mode into which particles 
or pairs can decay spontaneously. The respective decay 
rates Fi and r2 depend on the couplings set by the in- 
teraction constant, and on the correlations in the Bose 
liquids. To estimate Fi and r2 we resort to the univer- 
sal low-energy description of the Bose liquid in the filled 
mode in terms of a Luttinger liquid [37]. Within Fermi's 
golden rule, we find the decay rates to be [38]: 



T^iA.K) = U\lF2{K) 



1-1/2K* 



1/2-2/K' 



(3) 



Here Ai = A2 = AJ'^/A^^, relate to single- 

particle and pair tunneling between the upper and lower 
states, respectively, while A„ff = _)_ denotes 

the actual energy difference between these states. The 
functions Fi{K) w 47r and F2{K) « n/^/S are very 
weakly dependent on the Luttinger parameter K and 
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K* = K/{1 + Ssf' is an effective Luttinger parameter, 
renormalized by a non-universal phase shift due to the 
interaction between the quantum Hquid and the particle 
or pair transferred to the empty lower state (See ref. [39] 
and Supplementary Information). 

If we assume random independent decay events, then 
the density in the upper state obeys the equation p = 
— (Fi + 2T2)p^ = —Yp^- The decay is significant as long 
as the detuning between the tubes is less than the chemi- 
cal potential set by the interactions in the filled tube. At 
larger detuning, which exceeds the linear phonon spec- 
trum in that tube, the decay process is greatly suppressed 
due to phase space limitations. 

The density remaining in the upper state can now be 
found from the solution of the decay equation with the 
rates in Eq. (3) 



1 



2po 



dAr(A) 



(4) 



This result is valid in the limit of slow sweeps, when the 
time spent in the decay zone td = 2p/a is longer than 
the shortest response time of the Luttinger liquid ~ h/ p,. 
For shorter sweep times, the response is essentially that of 
free particles. Then the single-particle LZ formula Eq. (1) 
is expected to apply. 

The above arguments can be extended to the Mott in- 
sulator, where the opening of the gap in the excitation 
spectrum gradually eliminates the decay channels. This 
is taken into account by introducing the Mott gap Ami as 
a lower cutoff to the integral over V in Eq. (4). Deep in 
the Mott phase the decay vanishes completely as the sys- 
tem becomes a collection of effectively decoupled double 
wells. 

Taking all these points into account, we plot the ex- 
pected transfer efficiency tir as a function of the inverse 
sweep rate 27rJ^/|a| in Fig. 4b. The results have been 
calculated for different values of the axial 2;-lattice depth, 
taking into account the inhomogeneous filling of different 
tube pairs in the experiment (see Methods). Both exper- 
iment and theory show the same qualitative behavior. 

In addition to the breakdown of adiabaticity, we ob- 
serve an intriguing non-monotonous change in the trans- 
fer efficiency versus the axial z-lattice depth in the 
crossover between the slow and fast sweep regimes (0.5 < 
27rJ^/|a| < 1) (see Fig. 4a). The transfer efficiency nji 
first decreases and then increases with increasing z-lattice 
depth. A detailed investigation of this effect is plotted in 
Fig. 5a, where we show how n^j changes at a fixed value 
of 27rJ^/|a| with increasing z-lattice depth as the system 
evolves from a superfluid to a Mott insulator. The exper- 
imental data shows a pronounced minimum of nji prior 
to reaching the Mott transition. 

In order to explain this phenomenon, the loop struc- 
ture, caused by the non-linearity and plotted in Fig. 3b, 
needs to be taken into account. Let us consider the mean- 
field Hamiltonian (Eq. (2)) in the weak coupling regime, 
far from the Mott phase. Extremizing the mean-field en- 
ergy with respect to the inter-tube coherences leads to a 



non-linear version of the two state Schrodinger equation 
(see Supplementary Information). For a non-linearity 
77 > 1, the "energy levels" versus detuning develop a loop 
structure [42], as seen in Fig. 3b, which grows in width 
with increasing non-linearity. 

Such a loop structure leads to a dramatic dynamical 
evolution of the condensate once the resonance is crossed 
in the inverse sweep. First the condensate starts out far 
detuned on the right upper branch, which corresponds to 
a local extremum of the mean-field energy. In the course 
of a slow sweep the condensate moves along this branch 
as the detuning is decreased, and gets caught in the loop 
structure in a self trapped state. In a one-dimensional 
system such a state survives as a slowly decaying meta- 
stable state [40] . The edge of the loop is a spinodal point 
of the mean-field energy where the local extremal solu- 
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FIG. 5; Changing tlie loop size: influence of intra-tube corre- 
lations on transfer efficiency, a Measured transfer efficiency 
at a constant normalized sweep time 27rjVa = 0.53(3) as a 
function of the axial z-lattice depth and different inter-tube 
tunnel couplings. A pronounced minimum of the transfer ef- 
ficiency is observed prior to reaching the Mott transition in 
the quantum ladder. This minimum shifts towards higher z- 
lattice depths for decreasing inter-tube couplings J. b Phase 
diagram of the meta-stable upper condensate branch com- 
puted using mean-field theory (Eq. (2)). In the gray shaded 
area, the upper condensate branch exhibits a loop (Fig. 3b). 
The datapoints represent the minima in the measured trans- 
fer efficiency as obtained from a fourth-order polynomial fit to 
the data in a. They approximately coincide with the theoreti- 
cal prediction of the maximal loop size (maximal non-linearity 
in the mean-field picture). 
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tion disappears. When the detuning is decreased further, 
the condensate finds itself in a catastrophic scenario on 
a downward slope of the energy landscape, leading to a 
decay to the low-energy mode of the system via the dis- 
cussed emission of phonons into the ID quantum liquid 
(see Supplementary Information). 

In our case, the non-linearity 77 and therefore the size 
of the loop structure depends non-monotonously on the 
axial z-lattice depth. Far from the Mott transition the 
lattice acts to increase the effective on-site interaction 
and with it the effective non-linearity 77, which leads to 
growing loop size and enhanced probability of the con- 
densate to be self-trapped. As a consequence, nji first 
decreases with increasing z-lattice depth (see Fig. 5a). 
However, beyond a certain lattice depth the condensate 
$ begins to deplete significantly, causing a decrease of 77 
and hence an increase of nj^ until 7] vanishes close to the 
Mott transition where $ = 0. Fig. 5b charts the line of 
maximal non-linearity and the Mott phase boundary as 
a function of z-lattice depth and the inter-tube tunneling 
J. The former agrees well with the experimentally deter- 
mined position of the minima in transfer efficiency seen 
in Fig. 5a. Note that the decay mechanism still plays a 
role until deep in the Mott regime. Then, however, the 
phonon decay channel is shut off by the Mott gap and 
thus we recover the standard single-particle LZ behavior. 

Finally, let us comment on the probability for the sys- 
tem to become self-trapped in the loop structure for dif- 
ferent sweep rates. For very fast sweeps the behavior of 
the non-linear system is close to the linear one as the sys- 
tem is unlikely to be trapped in the loop structure. How- 
ever, for slower sweeps and larger non-linearities (leading 
to larger loop structures) the probability for self-trapping 
is strongly enhanced. We find that the crossover between 
the slow and fast sweep regimes, for strong non-linearity, 
takes place at inverse sweep rate J^/a* « 77"^/^ [41]. 



III. DISCUSSION AND CONCLUSIONS 

We have measured and analyzed novel dynamical phe- 
nomena in pairwise coupled ID Bose liquids, which con- 
stitute a direct generalization of the LZ dynamics in 
a many-body system. A particularly striking result is 
the sudden breakdown of adiabaticity seen in the inverse 
sweep, for which slower sweeps lead to lower transfer effi- 
ciency. We explained the results theoretically by treating 
the inter-tube coherence within a mean-field approxima- 
tion valid all the way to the strongly correlated regime. 
The analysis lead to an effective non-linear LZ problem 
similar to the one described in ref. [41, 42]. The break- 
down of adiabaticity in the inverse LZ sweep can then 
be attributed to self-trapping of the inter-tube coherence 
in a meta-stable solution, associated with a loop struc- 
ture in the energy levels, which emerges beyond a critical 
value of the non-linearity. We showed that the effective 
non-linearity decreases and eventually vanishes upon ap- 
proaching the transition to the Mott-insulating phase. A 



crucial new ingredient in the present setting is the cou- 
pling of the inter-tube coherence to the low energy exci- 
tations (phonons) of the ID Bose liquid. These fluctua- 
tions, treated within a quantum Luttinger liquid theory, 
provide the necessary dissipation mechanism, which al- 
lows for decay of the condensate from the meta-stable 
solution. 

Our work provides new insight into the quantum dy- 
namics of many-body systems far from equilibrium, and 
constitutes a step towards extending the LZ framework 
for treating such phenomena. An interesting scenario 
arises when our approach is applied to fermionic systems. 
The presence of an underlying Fermi surface can, in par- 
ticular, facilitate the observation of universal phenomena 
in the non-equilibrium quantum dynamics, such as pre- 
dicted for the Rabi oscillations in a similar setup con- 
sisting of a coupled tubes undergoing a sudden quench 
[38]. Another natural continuation of our work is to 
the dynamics of quantum spins in an optical lattice, for 
which the possibility of universal phenomena was sug- 
gested [43]. 



Methods 

Experimental cycle. We create the array of pair- 
wise coupled tubes by means of a bichromatic superlat- 
tice introduced in ref. [30]. The atoms are first loaded 
into a two-dimensional optical lattice with periodicity of 
As/2 = 765 nm along the a;-direction and Ay/2 = 420 nm 
along the y-direction. This is achieved by ramping up 
the lattices to Vs = AOE^ and Vy = 30 Ey, respec- 
tively, within 300 ms. In this way, we form a 2D array 
of tubes with about 100 particles in the central tubes. 
A short-period lattice along the x-direction with peri- 
odicity Aa;/2 = 382.5 nm allows us to split the x-lattice 
sites into double-wells. By ramping up this lattice to 
50£'if within 10 ms with the relative phase set to obtain 
highly tilted double- wells, we achieve loading of all atoms 
to left tubes alone. Any coupling in the y-direction is 
frozen out by ramping up the y-lattice to 50 within 
100 ms. During this ramp time, the z-lattice with wave- 
length A2 = 840 nm is ramped to the desired depth to 
tune the intra-tube tunneling Jy and the relative phase 
between x- and s-lattice is adjusted to the initial detun- 
ing A,; — El — Er, with El {En) the energy of the left 
(right) well {Ai/h = —5 kHz for the ground-state sweeps 
and Ai/h = 6.6kHz for the inverse sweeps). 

Before starting the LZ sweep, we set the inter-tube 
coupling J by rapidly ramping down the cc-lattice and 
thereby the double- well barrier in 200 fis. Then we per- 
form the LZ sweep by linearly changing the detuning be- 
tween left and right tube to a final value A f with a sweep 
time T — {Af ~ Ai)/a by changing the relative phase 
between the x- and s-lattice {Af/h — 5 kHz for ground 
state sweep and Af/h — —6.6kHz for inverse sweep). 
After the sweep is finished, we ramp up the barrier in 
200 /xs to quench the coupling again. By applying a site 
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resolving mapping method described in Refs. [30, 32], we 
measure the relative population in left and right tubes, 
respectively, in order to detect the transfer efficiency of 
the sweep. 

Bose-Hubbard ladder model. The double-tubes, 
with an optical lattice depth Vz > 5Ei^ along the tube 
axis, can be described in terms of a Bose-Hubbard ladder 
model: 

H = -J\\ E {blK + H.c) 

(ij),u=LM. 
i,v=L,R 

- E ^{n^R - n^L) + Jihlni^L + H-c)- (5) 

i 

By treating the inter-tube coherence in a mean-field ap- 
proximation we reduce from a two chain model to an 
effective one chain model for the quasi condensate mode 

al = V'L&L + ^Rh\R- 

= - J|| E + ^-c-) (6) 

i 

+ E [Ad^Lp - l^flP) - JiVL^R + c.c.)] a\a^. 

i 

This should be viewed as a variational Hamiltonian which 



defines a non-linear LZ problem for the inter-tube coher- 
ences V'LjV'-R 8.S described in the text. 

Transfer efficiency in inverse sweep. In order to 
compare this theory with the measurements we also need 
to account for the inhomogeneity in the experiment. Due 
to the transverse parabolic confinement, the experiment 
actually averages over tubes with different particle num- 
bers. The central tubes have high particle number, while 
more peripheral tubes have fewer particles. In quite a 
high fraction of the double-tubes, the chemical potential 
is lower than the level spacing between phonon states, 
and these act as zero dimensional double dots where the 
decay processes are absent. Taking into account those 
tubes we have the approximate formula for the slow 
sweep regime tir = (1 — vo)Pr + vq where vq is the frac- 
tion of the particles in the peripheral tubes. To test this 
hypothesis we repeated the experiment with varying par- 
ticle number. The expectation is that will decrease 
with increasing particle number N because more of the 
tubes will fill up considerably and enter the ID regime. 
In the experiment we indeed found that the saturation 
value of rifl, in the limit of slow sweeps, decreased with 
increasing N . 

We would like to thank Anatoli Polkovnikov, Giu- 
liano Orso, Christian Kasztelan and Ulrich SchoUwock 
for stimulating discussions. This work was supported 
by the DFG, the EU (STREP NAMEQUAM), DARPA 
(OLE program), AFOSR, DIP (EA and IB), and the ISF 
(EA). 
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S.I. 



Supplementary Information 

FEW-PARTICLE LANDAU-ZENER 
SWEEPS IN A DOUBLE DOT 



A simple picture for a multi-particle LZ sweep can be 
obtained when considering a double dot filled with N 
interacting particles. For strong repulsive interactions 
U ^ J, the ground state undergoes N well separated (in- 
dependent) tunneling resonances as the bias A is changed 
from —oo to oo [SI]. At each resonance, two Fock states 
\nL,nji) and \nL — l, ur + I) are coupled by single-particle 
tunneling. The corresponding coupling matrix element is 



J,^ = -J{nL,nR\a{dR\nL-l,nR + l) 
-\/nLinR + 1)J, 



(S.l) 



where the emerging prefactor > 1 stems from the bosonic 
statistics and + ur — N . We will now consider a lin- 
ear sweep with IS. — a-t starting from A(i = — oo) = — cxi 
and calculate the overall transfer efhciency from the left 
to the right well at t = oo. Herein, we will neglect the 
fc-th order tunnel couplings [k > 2) which occur between 
the excited states and are suppressed as [J/U)^^^ . The 
probability to stay adiabatic at a specific first order res- 
onance is given by substitution of Eq. (SI) into Eq. (1) 
as PLz{Jn„,(^). The probability to adiabatically follow 
the ground state throughout the whole sweep thus is 



N-l 



= n 



(S.2) 



1=0 



Taking into account the remaining right-well popula- 
tion at each diabatic crossing of a resonance, we obtain 
the transfer efficiency nR(o) of the sweep: 



N-l 



i=0 



LZ 



(S.3) 



In Fig. SI, we plot UR^a) for various numbers of par- 
ticles N in the double dot. The characteristic sweep rate 
ac increases gradually with N and we find an enhance- 
ment by a factor of 7 for = 10. 
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FIG. SI; Double-dot picture in the limit of strong interac- 
tions, a Plot of the eigenenergies versus A for a double dot 
filled with four particles and U/J = 30 The shaded regions 
mark the individual tunnel resonances crossed in the ground- 
state LZ sweep, b Plot of the right-well population after a 
LZ sweep in a double dot with strongly interacting particles 
{U 2> J) and rate a as calculated from Eq. (3). The trans- 
fer efficiency increases with increasing number of particles 
in the double dot. The inset shows the characteristic sweep 
rates with 71^(0^) = 1 — e^^. 



For smaller interactions U J , the tunnel resonances 
overlap and cannot be treated independently. We have 
performed numerical simulations of sweeps with U / J = 
0.5 using exact diagonalization. We plot nR{a) as ob- 
tained from these simulations for = 1 ... 20 in Fig. S2. 
The results show the increase of ac with N , where the 
enhancement for A^ = 10 is reduced to a factor of 2 as 
compared to the case of independent resonances. 

We note that the fidelity for an inverse sweep in a 
double dot is governed by the A^-th order coupling be- 
tween the states |A^, 0) and |0,A^) which is suppressed 
as (J/C/)^~^, yet finite. Therefore, it is always possi- 
ble to stay adiabatic for slow sweeps in the double dot 
system, even though the adiabaticity timescales might 
become unrealistically long. This is in contrast to the re- 
sults for pairwise coupled tubes presented in this article, 
where mean-field and Luttinger-liquid physics lead to a 
breakdown of adiabaticity for slow sweeps. 
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FIG. S2: Transfer efficiency for moderate interactions. We 
plot the simulated transfer efficiency n_R(a) for various num- 
bers particles and an interaction energy of i7 = 0.5J. Here, 
the tunnel resonances overlap. The enhancement of the char- 
acteristic sweep rate Uc (dashed line) with increasing A'^ per- 
sists, but becomes less strong compared to the case of inde- 
pendent tunnel resonances. 



S.II. ATOM NUMBER DISTRIBUTION IN THE 
OPTICAL LATTICE 

In our experiment, the presence of the magnetic trap 
and the Gaussian profile of the lattice beams leads to an 
overall harmonic confinement and correspondingly inho- 
mogeneous distribution of atom numbers in the tubes. 
The trap frequencies right after ramping up the initial 
2D optical lattice are Wj. = 27r x 45 Hz, ujy = 2tt x 40 Hz 
and = 27r X 59 Hz, respectively. Using the effective 
description of ID-gases (see e.g. [S2]) and an s-wave 
scattering length of as — 5.34 nm, we simulate the atom 
number distribution employing a local density approxi- 
mation. This results from the loading sequence described 
in the Methods section of the main article. For our typi- 
cal number of 9 x 10^ atoms, we find a maximum of 108 
particles per tube in the center (see Fig. S3) right after 
the initial rampup of the 2D lattice. 



S.III. THE BOSONIC LADDER 

Here we provide details for the theoretical descrip- 
tion of the coupled two-tube system with a focus on the 
Landau-Zener (LZ) passage discussed in the main arti- 
cle. The details given here cover two different aspects of 
the coupled Bose-liquids: (i) The energetics which can 
be captured in a mean-field description. The mean-field 
solution allows us to map the many-body system to an ef- 
fective two-state problem discussed by Wu et al. [S3, S4]. 
We show how to extract an adiabaticity criterion from 
this mapping, (ii) The dynamics of the collective excita- 
tions of the Bose liquids, which provide an effective bath 
for the LZ problem. In particular, we describe how to 
capture the decay processes responsible for breakdown of 
adiabaticity in the inverse sweep. 
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FIG. S3: Atom number distribution in the system, a Atom 
number distribution in the 2D lattice right after the initial 
loading. For our parameters and the total number of particles 
(9 x 10^), we get about 108 atoms in the central tube of the 
system and the occupation decreases with distance from the 
center. The mean value is 43 atoms per tube, b Histogram 
showing the total number of atoms residing in tubes with 
given occupation numbers. 



A. Non-linear Landau-Zener problem 

Before we analyze the Hubbard Hamiltonian on the 
two-leg ladder we review the physics of a simpler non- 
linear two-state model [S3, S4] described by the classical 
Hamiltonian: 

H[^L>fl] =-2JV'iVfl + A(|V;L|'-|Vfl|') (S.4) 
+ U{\i;L\^ + \M^). 

Here, A is the detuning which will be ramped as A = a-i. 
Note, that we consider [/ > 0. The two amplitudes tpR 
and ipL obey \i})r\^ + [■0l|^ = 1- 

For TpR^L real we parametrize V'-R — sin(7/2) and 
= cos(7/2). Fig. S4 shows the energy landscape as 
a function of the angle 7. For sufficiently strong non- 
linearity one finds four extrema of the energy near zero 
bias and only two further away. This is essentially the 
loop-structure discussed in [S3]. The loop structure ap- 
pears beyond a critical value of the non-linearity , i.e. 
r] = U/ J > 1. For r/ < 1 and A = 0, there are the two 
extrema corresponding to the bonding and anti-bonding 
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U/J=0 




FIG. S4: The energy of the non-linear Hamiltonian (S.4) as a 
function of the angle 7. The extrema mark the (meta-)stable 
solutions, a For no non-linearity {rj = 0), only two states 
are stable (dots), b Depending on the detuning either four 
(A — 0) or two (A = 1.5 J) states appear close to resonance. 

state of the non-interacting problem at 7 = tt/2 and 
7 = 37r/2, respectively. For rj > 1 the non-linearity is 
strong enough to give rise to four different extrema at 
zero detuning. When going away from A = 0, the detun- 
ing introduces a tilt and two of the extrema eventually 
merge and disappear. The resulting energy level struc- 
ture is given in Fig. 3b of the main article. 

The condensates 4'l,r obey the (Gross-Pitaevskii) dy- 
namics 

ihdtipL = d^iH{il^L,'^R), (S.5) 

prescribed by the Hamiltonian (S.4) together with the 
Poisson brackets {V'aiV'a} = 1- Based on this dynamics 
it is possible to compute a transition probability, using 
the adiabatic perturbation theory of classical mechanics 
[S5]. This yields [S4] F cx 1 - cxp(-ac/a) with ac = 
2TTj^q and 

q^- dx{l + x^)^/^[{l + x^)-^/^+r^f\ (S.6) 

^ Jo 



Note that g = 1 for 77 = 0, in accordance with the exact 
result. 



B. From a Bose-Hubbard ladder to two-mode 
dynamics 

In the main text we described a mapping from the 
two chain Hubbard model to the energetics of the non- 
linear two-mode LZ problem by employing a mean-field 
approximation for the inter-tube coherence. Here we will 
show that the dynamics is of the form (S.5) and derive 
the effective non-linearity 77 in this case. 

As described in the text, the key is to treat the inter- 
tube coherence in a single-mode approximation taking 
a uniform inter-tube wave-function for the condensate 

= V'L^i +'0-R^fl- The resulting Hamiltonian 

ff given 

by Eq. (6) in the Methods section is not exactly solvable. 
To make further progress we apply a mean-field (MF) 
approximation for both the axial and inter-tube modes. 
This of course misses the subtle (power-law) long range 
correlations in the one-dimensional Bose liquid, however 
it is not a bad approximation for the energetics. 

The MF theory can be easily formulated as a varia- 
tional approximation with j^*) taken to be a site factor- 
izable wave-function, which depends on the inter-tube 
coherences V'iiV'fl among other parameters. The time 
evolution, or the equations of motion for the variational 
parameters, are obtained from the time dependent vari- 
ational principle [S6, S7] 

C^,,^{^\ihdt-H\-^), (S.7) 

where H is the original ladder Hamiltonian (5). In the 
weak coupling limit U <C J\\n, the variational wave- 
function can be taken to be the coherent state 

J = e-l*l'+*('^-^I+'^«^I.)|vac), (S.8) 

which gives rise to the two-mode Gross-Pitaevskii dy- 
namics 

-J{rL^R + C.C.) + ^{\^L?-\^R\^) 

+ C/<I>^(|VL|' + |V'fl|')-2J||$2 (S.9) 

with non-linearity rj = C/|$p/J. The absence of dynam- 
ics in the axial degree of freedom is a result of constrain- 
ing the variational states to a uniform condensate $. In 
reality there are low-energy axial phonon excitations, the 
consequences of which will be discussed in the last sec- 
tion. 

For strong interactions we resort to the Gutzwiller MF 
theory, which consists in using the site-factorizable vari- 
ational states |^'0f,_H) ~ Hi fm\m)i. Here \m)i denotes a 
state with m Oj-bosons on site i. Very close to the Mott 
transition and in the Mott phase we may truncate to only 
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three occupation states, so that 



14- 



n 



cos(6'/2)|l) 



i(g/2)^ 
a/2 



(|0). + |2).) 



n cos(^?/2)[V'i6[ + Vfl6l,] 

i '- 



sm 



72 



I vac). 



We evaluate the effective Lagrangian C^a to leading order 
in small <!> = (^lail*) to obtain 

= + Vfl^lf) (S.IO) 

- J{rLi^R + c.c.) + AdV'Lp - IV^flH (S.ll) 
+ [/C$2(|V,i|4 + - (3 + 2\/2)J||$2 (S.12) 

+ Jc$2^ - |c|V^inv>,jpci>ci> + 4J||$V' (S.13) 



$4 



'J||(3 + 2y2) [/C(|^i|4 + |V^fl|4)^ 
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(S.14) 



The numerical factor C < 1 results from the truncation. 

In the lines (S.10)-(S.12) we recognize again the non 
linear two mode-dynamics for 'ipL,R with non-linearity 
rj = Ui^^ /J. In general this dynamics is coupled to the 
uniform dynamics of $ which correspond to a gapped 
condensate amplitude oscillation about the average value 
$ dictated by the theory [lines (S.12)-(S.14)] and de- 
pends on the distance from the Mott phase. For slow LZ 
dynamics we can therefore set $ and (p to their instanta- 
neous variationally determined values $ and — Q. 

To obtain the estimate for the critical ramp rate etc we 
compute the condensate fraction $ at A = for given 
U , J|| , and J . This gives the characteristic non-linearity 
parameter rj in the effective model. We then use rj to 
calculate the critical ac via (S.6). The result is shown in 
the insets of Fig. 2c and Fig. 2d in the main article. 



C. Phase diagram and maximal loop 

In this section we derive the phase diagram of the 
meta-stable Bose liquid in the upper state and deter- 
mine the point of maximal non-linearity for each given 
inter-tube coupling J as a function of the axial-lattice 
depth. This could be done, in principle, using the inter- 
tube MF approximation described above. However, to 
get the complete phase diagram including the region of 
very small J it is better to use a MF approximation which 
treats the inter-tube coupling exactly, keeping both the 
bonding and anti-bonding inter-tube states. This is done 
using the following family of site factorizable variational 



0)2 + O 2 o_ 



0.3 
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FIG. S5: Adiabatic phase diagram within the variational 
wave function approach (S.15). White denotes the Mott insu- 
lator, the color represents the expectation value of the kinetic 
energy. In the right panel we show in comparison the critical 
value of J|| obtained from the simpler mean-field approxima- 
tion described by Eq. (S.4) as a dashed line. 



states \'$) = Hi with 



|V'^) =cos0j cos||10),-fsin||01), \+ ^_ X 



1 00), 



sint( cos||20),+sin||02), )-|-cost|ll) 



(S.15) 



The Fock states \a/3)i denote occupation states at site i 
in the left (a) and right (/3) tube. This wave function, 
improves on the inter-tube MF, described in the previous 
section, by allowing for combinations of the states |a/3)j, 
which cannot be written using the operators Oj alone. 

The phase diagram computed using the states (S.15) 
is described as follows. First it is clear that for large de- 
tuning (A ^ J, U) one can neglect the lower tube (for 
the inverse sweep) and the phase diagram is that of a 
single chain. To understand the more interesting regime 
near resonance, we focus on the case A = for sim- 
plicity. The inter-tube hopping acts to reduce the effect 
of interactions, significantly shrinking the Mott phase. 
Fig. S5b displays a map of the superfluid order parameter 
I^Lp + \^R? in the space of J\\/U versus J/U. Fig. S5a 
shows the phase diagram for the case A = 0. These re- 
sults are in good qualitative agreement with Ref. [S8]. 
To get more precise results, Danshita et al. resorted to 
numerical DMRG calculations. 

The results shown in Fig. S5b are used in Fig. 5b. 
(i) The phase boundary shown in Fig. 5b is defined by 
the line where $ —J' 0. (ii) We took a contour of equal 
order parameter 4> = $o close to its maximal value as an 
estimate for the location of the maximal loop structure. 
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D. Relaxation due to axial phonons of the Bose 
liquid 

The crucial ingredient leading to the failure of adia- 
baticity in the inverse sweep, which was missing in the 
above MF considerations, is the existence of phonon 
modes in the ID quantum liquid. These modes facili- 
tate the decay of the inter-tube condensate wave- function 
from the meta-stable upper to the lower state via emis- 
sion of phonons. In this section we provide the details of 
the effective Hamiltonian used to theoretically account 
for this effect. 

We start by deriving the coupling of the two adiabatic 
modes. For the discussion of the quasi-condensate we 
only used the state = i/jl^Il + i^Rbln- To describe the 
coupling of the two adiabatic modes we also need the or- 



thogonal state &■ = — V'fi^Ii + V'L^I_R- As long as the loop 
structure does not develop, i.e. for weak non-linearity 
?7 <C 1, the two amplitudes 'ipL,R are approximately given 
by the solution of the single-particle double-well Hamil- 
tonian 




where A„ff = y/AJ-^ + A-^. This is valid for example in 
and close to the Mott phase. 

We can now express the interactions in terms of the 
operators and bi to obtain the coupling terms between 
the modes: 



We defined 



^ + ^b') + >'^2pa^Pb^ + M{par ~ Pbi) + H.c] + A2 [a\o^Ab^ + H.c] |. 



' -"t^ ^ X _i 8J2 , 2AJ 4J2 

Pai-a^a^, Ph, - b^b^, - i - ^12-^^, ^1-^2"' -^2-^- 

cl'f cff off cff 



(S.16) 



The coupling Hamiltonian (S.16) contains an intra-mode 
interaction oc Aq, an inter- mode interaction oc A12, and 
particle and pair transfer terms Ai, and A2, respectively. 

To understand the effect of the coupling on the 
Landau- Zener dynamics in the limit of slow sweeps, we 
make use of the effective long wavelength (bosonized) de- 
scription of the quantum liquid in the adiabatic upper 
mode. This is given by the Sine-Gordon model [S9]: 



27 



dx K\d,ipix)\^ + -\dJix)\'' 

A 



g dx cos 2-d{x) 



(S.17) 



Here Vg w I J\\n is the sound velocity with I the lattice 
constant. The coupling constant g and the Luttinger pa- 
rameter K are in general phenomenological parameters. 
In the superfluid phase K > 2 and it decreases with in- 
creasing lattice depth. The transition to a Mott insulator 
occurs at the universal value oi K = 2, below which the 
cosine term becomes relevant, giving rise to the Mott gap 
Ami. Note that the long wavelength theory is valid at en- 
ergies up to the UV-cutoff ujq ^ y/U J\\n above the upper 
inter-tube state. At the Mott transition ^ U . 

The dynamics of a particle, which decays to the lower 
inter-tube mode is described as a free particle at the bot- 



tom of an empty band 



■ dk n^k^ J 



(S.18) 



where m* « f? j J^t^ is the effective mass. 

The coupling Hamiltonian (S.16) can also be written 
in terms of the long-wavelength (bosonized) degrees of 
freedom. The most relevant term due to single particle 
tunneling maps to: 

ill = AiC/Aipi^/2 /■ dxe'^'^'''>Mx) + H.c. (S.19) 



Similarly for the two particle tunneling it is 

H2 - A2UX2Pa I dx e'^^^'''>Mx)Mx) + H.c. (S.20) 



Here, Ai and A2 are dimensionless non-universal con- 
stants of order 1. Finally, we must also consider the 
cross-interaction term A12, between the two modes. This 
gives rise to a "final-state interaction" between the quan- 
tum liquid in state a and the first particles transferred to 
the mode b. In this respect, the particles in the initially 
empty mode represent mobile scattering centers (impuri- 
ties) for the full mode. It is known that mobile impurities 
do not lead to back-scattering [SIO] at low energies. For 
this reason we can write the cross-interaction as 



H12 — Ai2UXi2Tr / dx pb{x)dx^{x 



(S.21) 
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Again A12 is a constant of order 1. The coupling (S.21) 
can be removed by a unitary transformation [Sll] 



U — exp 



iSs / dx' pi,{x')ip{x') 



(S.22) 



with ds — A12U X12K / hvsTT . This transformation leaves 
the Hamiltonians (S.17) and (S.18) invariant but trans- 
forms the operators as 



',(t) 



(t)/ 



?,(t), 



(S.23) 



We can now derive the decay rates Tii^2}{^t K) from 
the state a to 6 due to the tunneling terms (S.19) and 
(S.20) within a Fermi- golden-rule approach. For this we 
have to calculate the correlation functions 

Ci - {i^lix)^Ux)Mo)Mo)), 

C2 = (V'b^(x)Vil(a;)V'l(x)^l(x)V'b(0)v3fc(0)V'a(0)V'a(0)), 



for the Ai and A2 processes, respectively. Due to the 
"final-state- interaction" , the operators ipaix) acquire an 
extra exponent iSsip{x) which effectively changes the scal- 
ing dimension of the correlation functions Ci and C2. We 
can absorb this by rescaling the Luttinger parameter 



K* = 



K 



(S.24) 



Note that this rescaling is only relevant for the correla- 
tion functions Ci and C2 ■ For the description of processes 
in the isolated mode a, e.g., the discussion of instabilities 
of the Luttinger liquid, the relevant parameter is still the 
original K. The rates ri(2)(Aj-?^) given in the main ar- 
ticle are now obtained by standard manipulations [SI 2]. 
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